Shell-model study of the = 82 isotonic chain with a realistic effective hamiltonian 
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We have performed shell-model calculations for the even- and odd-mass A'" = 82 isotones, focus- 
ing attention on low-energy states. The single-particle energies and effective two-body interaction 
have been both determined within the framework of the time-dependent degenerate linked-diagram 
perturbation theory, starting from a low-momentum interaction derived from the CD-Bonn nucleon- 
nucleon potential. In this way, no phenomenological input enters our effective Hamiltonian, whose 
reliability is evidenced by the good agreement between theory and experiment. 

PACS numbers: 21.60.Cs, 21.30.Fe, 27.60.-t-j 



I. INTRODUCTION 

The N = 82 region has long been the subject of shell- 
model studies owing to the strong doubly magic char- 
acter of ^■^^Sn. In this context, several calculations have 
been performed employing purely phenomenological two- 
body effective interactions Clearly, a main issue in 
the study of this region is that it provides the oppor- 
tunity to investigate the effect of adding protons to a 
doubly-magic core over a large number of nuclei (from 
mass number A=133 to A=152). Comprehensive stud- 
ies of this kind were conducted some twenty years ago 
m Refs. 0] and In the former a simple pairing 
Hamiltonian and empirical single-particle (SP) energies 
were used, while in the latter the two-body matrix el- 
ements and single-particle energies were determined si- 
multaneously by a least squares fit to the ground-state 
and excited-state energies drawn from all experimentally 
studied = 82 nuclei. 

Starting in the mid 1990s, however, shell-model cal- 
culations employing effective interactions derived from 
realistic nucleon-nucleon (NN) potentials have been per- 
formed in the N = 82 region, which have generally 
yielded very good results iJ, H, 0, H] ■ Some of these 
studies d, 13, 13] have focused attention on the evolution 
of the low-energy properties of the A^ = 82 isotones as a 
function of A, showing that the observed behavior is well 
reproduced by the theory. In Refs. [1,01 the shell-model 
results were also compared with those of a quasiparti- 
cle random-phase approximation calculation leading to 
the conclusion that the low-lying states are equally well 
described by these two approaches. 

A main merit of the above mentioned realistic shell- 
model calculations is, of course, that no adjustable pa- 
rameter appears in the matrix elements of the two-body 
effective interaction. In all of them, however, the single- 
particle energies have been taken from experiment. Actu- 
ally, the calculational techniques used in these studies are 
based on the time-dependent degenerate linked-diagram 
perturbation theory [9'| , which provides an effective shell- 



model hamiltonian, i?eff, containing both one- and two- 
body components. Usually, however, the one-body com- 
ponents, which represent the theoretical SP energies, are 
subtracted from iJcff and replaced by those obtained from 
the experimental spectra of nuclei with one valence nu- 
cleon [10l| . This procedure has indeed led to a very good 
description of nuclear properties in different mass regions 

m- 

Recently, we have performed fully realistic shell-model 
calculations employing both theoretical SP energies and 
two-body interactions in studies of p— and sd— shell nu- 
clei [l3, 113 ■ these studies we have renormalized the 
high-momentum repulsive components of the bare NN 
potential Vnn by way of the so-called Mow-k approach 
m, EBl, which provides a smooth potential preserving 
exactly the onshell properties of the original Vn n up to 
a cutoff A. The effective hamiltonian has then been de- 
rived within the framework of the time-dependent degen- 
erate linked-diagram perturbation theory. In both papers 
(T^fisj a cutoff momentum A = 2.1 fm"'^ was employed. 
This corresponds to the laboratory energy Eiab — 350 
MeV, which is the inelastic threshold of the A^A^ scat- 
tering. The results of these studies compare well with 
experiment and other shell-model calculations [H, [l3| . 

On these grounds, we have found it interesting to carry 
out a similar shell-model study of the A^ = 82 chain, fo- 
cusing attention on the low-lying states of both even- and 
odd-mass isotones. In this case, the SP energies and two- 
body effective interaction are determined from a Viow-k 
derived from the high-precision CD-Bonn potential [l^ 
with a cutoff momentum A = 2.6 fm~^. This value of A, 
which is somewhat larger than that used in our studies of 
the light nuclei, is needed to obtain a reasonable descrip- 
tion of the experimental SP spectrum of the medium- 
heavy nucleus ^^'^Sb. This is mainly due to the fact that, 
at variance with the model spaces for light nuclei, in this 
case the 50-82 major shell contains the intruder 0/iii/2 
state, whose relative energy turns out to be strongly sen- 
sitive to the value of A. More precisely, for small values 
of A this state lies far away from the other SP levels of 



the 50-82 shell, while the correct structure of this shell is 
restored when increasing the value of A f]3\ . 

The dependence of the results on A would obviously 
vanish when complementing the Viow-k effective two- 
body interaction with different three- and higher-body 
components for each value of A. However, it is at present 
a very hard task, from the computational point of view, 
to take into account two- and three-body forces on an 
equal footing in a third-order perturbative calculation for 
medium- and heavy-mass nuclei. Therefore, the value of 
the cutoff A may be considered as a parameter that we 
have chosen so as to reproduce satisfactorily the exper- 
imental SP spectrum. For the sake of completeness, it 
should be mentioned that the trend of our calculated SP 
energies suggests that a larger value of the cutoff could 
even improve the agreement with experiment. However, 
this would amount to include higher-momentum compo- 
nents, which would rapidly deteriorate the convergence 
properties of the perturbative expansion used to derive 
the effective shell-model hamiltonian. In the following 
section we show that our choice of A is a reasonable com- 
promise to assure both good perturbative properties and 
quality of the results. 

The paper is organized as follows. In Sec. II we give 
a brief outline of our calculations, focusing attention on 
the convergence properties of the effective hamiltonian. 
Sec. HI is devoted to the presentation and discussion of 
our results, while some concluding remarks are given in 
Sec. IV. 



II. OUTLINE OF CALCULATIONS 

Our goal is to derive an effective hamiltonian for shell- 
model calculations in the proton 2sld0g0h shell, which is 
the standard model space to describe the spectroscopic 
properties of TV = 82 isotones. Within the framework 
of the shell model, an auxiliary one-body potential U is 
introduced in order to break up the hamiltonian for a 
system of A nucleons as the sum of a one-body term Hq, 
which describes the independent motion of the nucleons, 
and a residual interaction Hi: 

H = + y v.f'' = T + v = 

= {T + U) + {V~U) = Ho + Hi . (1) 

Once Hq has been introduced, a reduced model space 
is defined in terms of a finite subset of Hq^s eigenvectors. 
The diagonalization of the many-body hamiltonian ([1]) in 
an infinite Hilbert space, that it is obviously unfeasible, 
is then reduced to the solution of an eigenvalue problem 
for an effective hamiltonian H^s in a finite space. 

In this paper, we derive i?eff by way of the time- 
dependent perturbation theory Namely, HcS is 
expressed through the Kuo-Lee-Ratcliff (KLR) folded- 
diagram expansion in terms of the vertex function Q-box, 



which is composed of irreducible valence-linked diagrams 
[20I [2H . We include in the Q-box one- and two-body 
Goldstone diagrams through third order in Hi. The 
folded-diagram series is summed up to all orders using 
the Lee-Suzuki iteration method ^22i |. 

The hamiltonian Hcs contains one- body contributions, 
whose collection is the so-called S-hox . As mentioned 
in the Introduction, in realistic shell-model calculations 
it is customary to use a subtraction procedure, so that 
only the two-body terms of -ffeff, which make up the ef- 
fective interaction VeS, are retained while the SP ener- 
gies are taken from experiment. In this work, we have 
adopted a different approach employing SP energies ob- 
tained from the iS-box calculation. In this regard, it is 
worth pointing out that, owing to the presence of the 
—U term in Hi, [/-insertion diagrams arise in the Q- 
box. In our calculation we use the harmonic oscillator 
(HO) potential, U — imw^r^, and take into account all 
[/-insertion diagrams up to third order. The oscillator 
parameter is ftw = 7.88 MeV, according to the expres- 
sion m = - 25^-2/3 foj. ^ ^ 1^2. 

Let us now outfine the Viow-k approach [13, [3 to the 
renormalization of Vnn- The repulsive core contained in 
Vat AT is smoothed by integrating out the high-momentum 
modes of Vnn down to a certain cutoff A. This integra- 
tion is carried out with the requirement that the deuteron 
binding energy and phase shifts of Vnn are preserved 
by Viow-k; which is achieved by the following T-matrix 
equivalence approach. We start from the half-on-shell T 
matrix for Vjvat 

T{k',k,k^) = VNN{k',k) + 

V q^dqVNN{k',q)^^-^T{q,k,k^) , 

where V denotes the principal value and k, k' , and 
q stand for the relative momenta. The effective low- 
momentum T matrix is then defined by 

Tiow-kip' ,P,P^) = Viow-kip',p) + 

V / q^dqViow-kip ,q)^ ^71ow-k((7,P,P^) , 

Jo P - q 

where the intermediate state momentum q is integrated 
from to the momentum space cutoff A and {p' ,p) < A. 
The above T matrices are required to satisfy the condi- 
tion 

T{p',p,p^) = Tiow^k{p',P,p'^) ; {p',p) < A. (4) 

The above equations define the effective low- 
momentum interaction Viow-k, and it has been shown 
[isl l that their solution is provided by the KLR folded- 
diagram expansion [23| mentioned before. In addition 
to the preservation of the half-on-shell T matrix, which 
implies preservation of the phase shifts, this Viow-k pre- 
serves the deuteron binding energy, since eigenvalues are 
preserved by the KLR effective interaction. For any value 
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of A, the Viow-k can be calculated very accurately us- 
ing iteration methods. Our calculation of Vinw-k is per- 
formed by employing the method proposed in [2J| , which 
is based on the Lee-Suzuki similarity transformation [2^ . 




FIG. 1: Diagonal r = 0+ TBME of V^s as a function of 
A'^max (see text for details). 



As mentioned in the Introduction, our Mow-k has been 
derived from the CD-Bonn NN potential with a cutoff 
momentum A = 2.6 fm~^, and for protons the Coulomb 
force has been explicitly added to Viow-k- 

A discussion of the convergence properties of H^s is 
now in order. For the sake of clarity, we first consider 
the two-body matrix elements (TBME) of Vcs and then 
the SP energies. 



TABLE L Q-box J'" = 0^ two-body matrix elements at third 
order in Viow-k (in MeV) compared with those obtained by 
calculating the Fade approximant [2|1]. 



Configuration 


3rd order 


Fade [2|1] 


((057/2)1Vq5„,|(057/2)"> 
{{Qg7/2f\VQbo.\{ld^/2f) 
{(Off7/2)'|VQi,o.|(ld3/2)') 
((057/2)'|Vgi,o.|(2Sl/2)2> 


-0.8624 


-0.8757 


-1.0022 


-1.0196 


-0.6948 


-0.6999 


-0.4642 


-0.4662 


((0<77/2)'|l^Q6o.|(0feii/2)') 
((ld5/2)'iVQi,..|(ld5/2)'> 
((ld5/2)'[VQ6o.|(ld3/2)') 
m^/2f\VQbo.\{2s^/2?) 


2.0532 


2.058 


-0.7069 


-0.7138 


-1.7611 


-1.8191 


-0.6807 


-0.7555 


((ld5/2)'|VQi,o.|(0/ln/2)') 
((ld3/2)'lVQi,..|(ld3/2)') 
((ld3/2)'|l^Qi,o.|(2Sl/2)') 

((ld3/2)'|V^Qi,o.|(0/lll/2)') 
((2Si/2)'|VQi,„,|(2Si/2)') 

((2Si/2)'|Vgi,o.|(0/lll/2)') 
((0/lll/2)'l Vote. 1(0/111/2)') 


1.2028 


1.2081 


0.0838 


0.0447 


-0.4817 


-0.4857 


1.0457 


1.0459 


-0.882 


-0.4686 


0.658 


0.6586 


-0.8025 


-0.8391 



In Fig. m we report the diagonal J"" = 0+ TBME 
of T4ff as a function of the maximum allowed excitation 
energy of the intermediate states expressed in terms of 



the oscillator quanta A^max- We have chosen the J'^ = 0+ 
matrix elements because they are the largest ones in the 
entire set of VeS TBME and play a key role in determining 
the relative spectra and ground-state (g.s.) properties 
of the even A^ = 82 isotones. From Fig. [T] it is clear 
that our results have practically achieved convergence at 

A^max — 15. 

As regards the order-by-order convergence, an estimate 
of the value to which the perturbative series should con- 
verge may be obtained by using Fade approximants. In 
Table HI we compare all the Q-box J'^ = 0+ TBME caL 
culated at third order with those given by the [2|1] Fade 
approximant [25l | 



^Qbox 



^Qboxl ^Qbox 



(5) 



^Qbox being the 7ith-order contribution to the J'^ = 0+ 
TBAIE in the linked-diagram expansion. From Table U 
we see that our third-order results are in good agreement 
with those from the [2|1] Fade approximant, the differ- 
ences being all in the order of few tens of keV with the 
exception of the diagonal matrix element for the (2si/2)^ 
configuration. This indicates a weak dependence of our 
results on higher-order Q-box perturbative terms. 




FIG. 2: (Color online) Theoretical relative SF energies as 
a function of Amax (see text for details). The experimental 
spectrum of ^^^Sb is also reported. 



In Fig. [2] we report the relative SP energies calculated 
at third order in Hi as a function of A^max- In the same 
figure, the results obtained by calculating the Fade ap- 
proximant [2|1] are also reported (red lines). Once again, 
the agreement between third-order results and those from 
this approximant points to a good perturbative behavior 



of the calculated relative SP energies. We also see that 
an A'niax — 15 calculation yields satisfactorily converged 
values for the relative SP energies. 

A different situation occurs for the absolute SP ener- 
gies. Actually, the calculated energy of the 1^5/2 level 
relative to ^^zsn is -4.468, -5.225, and -5.954 MeV for 
A^max=ll, 13, and 15, respectively. This shows that 
there is no sign of convergence of the absolute SP en- 
ergies as a function of the number of intermediate states. 
Moreover, the Pade approximant [2|1] of the 1^5/2 en- 
ergy for A'max = 15 is equal to -7.262, showing that also 
the order-by-order convergence of the absolute SP ener- 
gies is unsatisfactory. This poor perturbative behavior 
may be traced to the chosen value of the cutoff momen- 
tum, A = 2.6 fm~^, which is substantially larger than 
the standard one, A ~ 2.1 fm~^. However, the inac- 
curacy of the absolute SP energies affects only the g.s. 
energies of the nuclei considered. In fact, the calculated 
spectroscopic properties (energy spectra, electromagnetic 
transition rates) are certainly reliable, based on the good 
convergence of the relative SP energies and TBME. 
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FIG. 3: Experimental and calculated ^'^''Te spectra. 



III. RESULTS AND DISCUSSION 

We have per formed calculations, using the Oslo shell- 
model code [26| , for the even-mass N = 82 isotones up to 
^^''Hf, which is the last known nucleus belonging to this 
chain, and for the odd-mass ones up to ^''^Ho, the last 
odd-mass nucleus with a well-established experimental 
low-lying energy spectrum. This is the longest isotonic 
chain approaching the proton drip line, and therefore pro- 
vides an interesting laboratory to study the evolution of 
nuclear structure when adding pairs of identical particles. 



TABLE II: Calculated SP relative energies (in MeV), com- 
pared with the experimental spectrum of ^'''^Sb [JJ. The val- 
ues in parenthesis are the absolute SP energies with respect 
to the doubly closed ^''^Sn. 



nlj 



Calc. 



Expt. 



O57/2 



Id 
Id 



5/2 

'3/2 



2Sl/2 
0/lll/2 



0.000 (-6.999) 
1.045 
2.200 
2.006 
3.645 



0.000 (-9.663) 
0.962 
2.440 

2.793 



In Table |TT] our calculated SP energies are reported 
and compared with the experimental spectrum of ^^'^Sb. 
We see that the latter is on the whole reasonably well 
reproduced by the theory. However, while the calculated 
positions of the 1^5/2 and 1^3/2 levels come very close 
to the experimental ones, the energy of the O/111/2 level 
is overestimated by about 0.85 MeV. As we shall discuss 
later in this section, the energy of the 2si/2 level, for 
which there is no experimental information, appears to 
be somewhat underestimated. 



A strong test for our effective hamiltonian is given by 
the calculation of the energy spectrum of ^^"^Te, since the 
theory of the effective interaction is tailored for systems 
with two valence nucleons. From Fig. [31 where the ex- 
perimental [13] and calculated ^'^^Te spectra are reported 
up to 3.5 MeV excitation energy, we see that a very good 
agreement is indeed obtained. 




10 12 14 16 18 20 22 



FIG. 4: (Color online) Experimental and calculated ground- 
state energies per valence proton for A' = 82 isotones from 
A — 134 to 154. Zvai is the number of valence protons. 



In Fig. m we show the calculated and experimental 
[2^ g.s. energies (relative to the ^^^Sn core) per valence 
proton as a function of the number of valence particles 
Zvai of even-mass isotopes. We see that the experimental 
and theoretical curves are practically straight lines hav- 
ing the same slope, while being about 2.4 MeV apart. 
This discrepancy is essentially the same as that existing 
between the theoretical and experimental g.s. energies of 
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^^■^Sb (see Table |TT|. This confirms the reliability of our 
SP spacings and TBME, since the pattern of the theo- 
retical curve depends only on these quantities. 



5- 2.0 




for the J'^ = 2+ and 4+ states, for which the discrepan- 
cies do not exceed 350 keV. However, it should be noted 
that according to our calculations the proton subshell 
closure, which experimentally occurs at "^Gd, is not re- 
produced. As mentioned before, this can be traced to the 
theoretical position of the 2si/2 orbital, which appears to 
be too low in energy. This is confirmed by the fact [3| 
that this level has to be placed at 2.8 MeV in order to re- 
produce the experimental energy of the J'" = at 2.150 

137 



MeV in '''Cs, which is predominantly of SP nature [29(1 . 
We have verified that, if the 2si/2 SP level is placed at 2.8 
MeV, we obtain that the 2^ excitation energy in ^"'^Gd 
raises from 1.716 MeV to 1.976 MeV, while in ^^^Dy it 
decreases from 1.827 to 1.758 MeV, thus reproducing the 
observed subshell closure. 



FIG. 5: (Color online) Experimental and calculated excitation 
energies of the yrast J'^ = 2"*" states for N — 82 isotones. 
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FIG. 6: (Color online) Experimental and calculated excitation 
energies of the yrast = 4'^ states for N — 82 isotones. 




FIG. 7: (Color online) Experimental and calculated excitation 
energies of the yrast J'^ = 6'^ states for N = 82 isotones. 



In Figs. [5l[6l and[7]we report, as a function of the mass 
number A, the experimental and calculated excitation 
energies of the 2+, 4+, and 6+ yrast states, respectively. 
The experimental behavior is well reproduced, expecially 



TABLE III: Experimental and calculated B{E2;2t Of). 
The reduced transition probabilities are expressed in W.u.. 



Nucleus 



Calc. 



*Te 



"«Xe 
i^^Ba 

"^Nd 
"*Sm 



Expt. 



5.5 
9.06 
11.1 
14.5 
16.26 
16.6 



6.3 ±2.0 
16.6 ± 2.4 
10.8 ±0.5 

13.8 ±0.3 
12.03 ± 0.22 

11.9 ± 0.4 



To have a more complete test of the theory, we have 
also calculated the B{E2;2^ — > 0^) transition rates up 



to 



144 



Sm employing an effective operator obtained at 
third order in perturbation theory, consistently with the 
derivation of i/cft- Our results are reported and com- 
pared with the experimental data in Table IIIII We see 
that the agreement is quite good, providing evidence for 
the reliability of our calculated effective operator which 
takes into account microscopically core-polarization ef- 
fects. 

Let us now come to the odd-mass isotones. In Figs. [51 
l^fTUl andfTDwe report, as a function of the mass number 
A, the experimental and calculated energies of the 

, 2 J and — yrast states relative to the ^ yrast 
state, respectively. The experimental behavior is well re- 
produced for the 1^, I"*", and ^ states. For the latter, 
however, most of the calculated energies are higher than 
the experimental ones by about 300-400 keV, which re- 
flects the theoretical overestimation of the value of the 
0^*11/2 SP energy. Less satisfactory is the comparison 
between the calculated and experimental behavior of the 
^'^ yrast states, which is a further confirmation that the 
theoretical 2si/2 level in ^^'^Sb lies too low in energy. As 
a matter of fact, by performing a calculation within the 
seniority scheme up to d = 3, we have found that in ^^^Cs 



and ^■^^La the lowest 



1 + 



2 states, which are quite well re- 
produced by the theory, have significant seniority w = 3 
components. From A = 141 on, the experimental (5^)1 
levels are predominantly of SP nature, as testified by the 
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spectroscopic factors [27j, and the corresponding calcu- 
lated states, which are dominated by v = 1 components, 
underestimate their energies. 




133 135 137 139 141 143 145 147 149 
A 



FIG. 8: (Color online) Experimental and calculated excitation 
energies of the yrast J'^ = | states relative to J'^ = (| )i 
for odd-mass N = 82 isotones. 
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FIG. 9: (Color online) Experimental and calculated excitation 
energies of the yrast J'^ = | states relative to J'^ = (| )i 
for odd-mass N = 82 isotones. 




133 135 137 139 141 143 145 147 149 
A 



FIG. 10: (Color online) Experimental and calculated ex- 
citation energies of the yrast — ^ states relative to 
J" = (f ^)i for odd-mass N = 82 isotones. 
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FIG. 11: (Color online) Experimental and calculated exci- 
tation energies of the yrast J'^ = ^ states relative to 
— (|^)i for odd-mass TV = 82 isotones. 

As pointed out in the Introduction, the main purpose 
of this study has been to test the use of a fully realistic 
shell-model hamiltonian for the description of medium- 
heavy mass nuclei. In this context, it should be noted 
that, aside from adopting experimental single-particle en- 
ergies, earlier calculations on the N=82 isotones differ 
from the present ones in other respects, like the starting 
NN potential and the renormalization procedure. For 
instance, while we employ the CD-Bonn potential renor- 
malized through the V|ow-k procedure, in the works of 
Refs. d, Q use is made of the Bonn A potential and the 
Brueckner G-matrix approach. A comparison between 
the results of earlier works and the present one would 
not therefore be very meaningful. It is worth mention- 
ing, however, that the agreement between experiment 
and theory achieved in this paper is on the whole com- 
parable with that obtained in the previous calculations 
employing realistic effective interactions. 



IV. CONCLUDING REMARKS 

As discussed in detail in the Introduction, the N=82 
isotonic chain has long been considered a benchmark for 
shell-model calculations. This has resulted in a number 
of theoretical works which have practically all led to a 
good description of low-energy properties of these nuclei, 
evidencing the strong doubly magic character of ^'^^Sn. 

A main step towards a microscopic shell-model descrip- 
tion of nuclear structure has been the use of two-body ef- 
fective interactions derived from the free nucleon-nucleon 
potential. This approach, in which no adjustable pa- 
rameter is involved in the calculation of the TBME, has 
been successfully applied to the N=82 isotones in the 
last decade, showing the ability of realistic effective in- 
teractions to provide an accurate description of nuclear 
structure properties. 

However, some important theoretical questions still re- 
main open. In particular, these concern the calculation 
of the single-particle energies and the role of three-body 
correlations. In the present paper, we have tried to inves- 
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tigate the former issue by constructing a realistic effective 
shell-model hamiltonian for the = 82 isotones, where 
both the SP energies and TBME have been obtained 
starting from a Viow-k derived from the CD-Bonn poten- 
tial. This is a natural extension of our studies of shell 
nuclei [12] and oxygen isotopes [13] ■ Here, we have shown 
that our effective hamiltonian yields results which are in 
quite good agreement with experiment along the whole 
chain of the = 82 isotones. As regards the calculated 
single-particle energies, significant discrepancies with the 
experimental values occur only for the O/111/2 and 25^^/2 
levels. Note that for the latter no direct comparison is 
possible since the single-particle 2si/2 state is still miss- 
ing in the experimental spectrum of ^■^'^Sb. However, 
we have verified that increasing the calculated value by 
0.8 MeV the proton subshell closure occurs at ^^^Gd, in 



agreement with the experimental findings. 

At this point, the question naturally arises as to how 
contributions from three-body forces would modify the 
present results, in particular as regards the single-particle 
energies. While this remains a major task for future in- 
vestigations, we feel that the present study paves the way 
toward a better understanding of the microscopic foun- 
dations of the nuclear shell model. 
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